Screening of hydrodynamic interactions for polyelectrolytes in salt solution 
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We provide numerical evidence that hydrodynamic interactions are screened for charged poly- 
mers in salt solution on time scales below the Zimm time. At very short times, a crossover to 
hydrodynamic behavior is observed. Our conclusions are drawn from extensive coarse-grained com- 
puter simulations of polyelectrolytes in explicit solvent and explicit salt, and discussed in terms of 
analytical arguments based on the Debye-Hiickel approximation. 
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Macromolecules of biological relevance such as DNA 
are often highly charged polyelectrolytes. Under physi- 
ological conditions, these molecules are dissolved in salt 
solutions with Debye lengths of less than 1 nm. From 
the point of view of statics, the electrostatic interac- 
tions are thus screened: Large DNA strands in physi- 
ological buffers have similar static properties than reg- 
ular self-avoiding chains with short range interactions. 
The dynamical properties are more complex. For neu- 
tral chains, solvent-mediated hydrodynamic interactions 
influence the mobility and the internal modes of the 
molecules. The signature of these hydrodynamic inter- 
actions is dynamical "Zimm scaling"—. For example, the 
mobility fi and the diffusion constant _D of a self-avoiding 
Zimm chain scale as ^TV^" with the chain length 

N, where = 0.588 is the well-known Flory exponent. In 
contrast, the same quantities scale as in a "Rouse" 
chain, where hydrodynamic interactions are absent or 
screened. 

For polyelectrolytes like DNA, the interaction with 
the counterions and salt ions complicates the situation. 
The polyelectrolyte chain is surrounded by an oppo- 
sitely charged ion cloud, which drags behind and gen- 
erates additional friction^. This "relaxation effect" re- 
duces the mobility, but the dominating behavior for large 
chain lengths is still predicted to be Zimm-like, and in- 
deed, the scaling exponents reported in experiments (e.g., 
D - iV-°-^^ for double stranded DNA in Ref. i) are 
closer to the Zimm than to the Rouse exponent. 

On the other hand, it has long been known that 
the electrophoretic mobility in an applied electrical field 
scales like /ig ~ for long chains, i.e., it exhibits 

Rouse-like behavior—. The physical explanation for this 
effect is rather intuitive: The electrical field acts both on 
the polyelectrolyte and the surrounding ion cloud. The 
net momentum transferred to the solvent by all parti- 
cles is thus zero, and the hydrodynamic interactions are 
screened as a result. One can derive a screened Oseen 
tensor, which describes hydrodynamic interactions be- 
tween monomers in a homogeneous electrical field-j^i. 

Hence it appears that electrostatic and hydrodynamic 
screening are closely related to each other, and that elec- 
trostatic screening may entail hydrodynamic screening if 
the dynamics is driven by electrical fields. Electrostatic 
forces drive electrophoresis, but also to some extent the 



internal motion of polyelectrolytes, since they dominate 
the nonbonded interactions between monomers. The 
question thus arises whether hydrodynamic screening can 
also be observed in the absence of external fields. 

This question is addressed in the present letter. We 
consider a polyelectrolyte in salt solution, without exter- 
nal fields, and study the dynamics of internal motions 
in the chain. The dynamical scaling - Zimm vs. Rouse 
- gives us information on the influence of hydrodynamic 
interactions. We use a generic coarse-grained model both 
for the polyelectrolyte and the solvent. Recently, it has 
been shown that such models can account quantitatively 
for the dynamics of real polyelectrolytes in solutionSiiS. 

On general grounds, we already know that the dynam- 
ics on very large time scales, where chain diffusion dom- 
inates, must be Zimm-like: By virtue of the Einstein re- 
lation, the diffusion constant is related to the mobility 
of the chain in a sedimentation field, which acts only on 
the monomers and not on the ions. The interaction be- 
tween monomers and ions leads to a relaxation effect, 
as discussed above, but not to hydrodynamic screen- 
ing. Diffusion-dominated behavior sets in at the "Zimm 
time", i.e., the time scale on which single monomers fol- 
low the center of the mass motion of the whole chain. 
In the present work, we thus focus on time scales below 
the Zimm time, where the dynamics is governed by the 
internal modes of the chain. 

We study the dynamics of single bead-spring chains 
immersed in a solution of solvent particles and ions (salt 
ions and counterions) by dissipative particle dynamics 
(DPD) simulations^. Ions and monomers are hard-core 
particles with diameter cr, which interact with the po- 
tential Kc = 4e((r/cr)i2 - (r/cr)S) at r < cr. In ad- 
dition, the ions and a fraction of the monomers carry 
single-valued charges ±e and interact via Coulomb po- 
tentials. The monomers in the chain are connected by 
finitely extensible nonlinear elastic (FENE) springs with 
the bond potential Ko„d = 1/2 k R'^ln{l - {r/Rf ). 
The solvent particles provide the fiuid background and 
have no conservative interactions. All particles inter- 
act with dissipative DPD forces, which have the usual 
formi FfJ?^-^ = YP, + Yp, with a viscous contribution 
F,^ = — 7ti;(r. )(f.. • v .)f.. and a random force = 
y^2jk^ Tuj {r.. )r.. Xij ■ Here v is the velocity difference 
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FIG. 1: Static structure factor S{k) of the charged chain in 
salt solution (solid line) compared to that of an uncharged 
chain (dashed line). The inset shows the mean square dis- 
placement of the central monomer relative to the center of 
mass of the chain ((Ar)^) as a function of the rescalce time 
Dt with the diffusion constant D. 



and r.^ the distance between the two particles i and j, 
= ^{r) = 1 - {r/a) for r < a, and x,, is 

a Gaussian distributed white noise with mean zero and 
variance one. All particles have the same mass m. The 
natural units in our system are thus the length a, the 
thermal energy k^T, and the time unit r = ^mjkgTa. 

Specifically, we consider chains of length = 50 with 
charges on every second monomer, in solutions with ion 
concentration p. — O.lcr"^. The other model parameters 
are e = ak^T, e = fc^T, R = 1.5a, k = 25k^T/a^, 
and 7 = bk^Tr/a"^, and the density of solvent parti- 
cles is = 3a~^. With these parameters, the motion 
of the ions in solution is diffusive after an initial ballis- 
tic regime of length ~ lOr. The Debye length in the 
solution is A^^ = 0.89<t, the kinematic viscosity of the 
fluid is 7] ~ \.2AkgTT / , and the ions have the mobil- 
ity C, = OAAa^ /rk^T. The mean radius of uncharged 
chains of length TV = 50 is i?^ „ = A.TAa, and they diffuse 
with the diffusion constant = Q.021a^ /t. In compar- 
ison, the charged chains are just slightly more swollen, 
Rg ^ = 5.35cr, but their diffusion constant is significantly 
reduced due to the relaxation effect, Dc — 0.013(t^/t. 
The simulations were carried out in cubic simulation 
boxes with system size L = 25a and periodic bound- 
ary conditions, using the time step At = O.Olr. The run 
lengths were 4 million time steps after an equilibration 
time of 2 million steps. 

Fig. [T] compares the static structure factor for such 
a charged chain in salt solution to that of an uncharged 
chain with the same length. The charged chain is slightly 
stretched, otherwise the two structure factors are very 
similar. Algebraic scaling behavior S{k) oc fc^" is ob- 
served in the wavevector range 1/R < k < l/ftg, where 
flg 1.7a is the distance between charged monomers in 
the chain. The value of v is found to be = 0.67 for 
the uncharged chain, and v = 0.7 for the charged chain. 
These values are larger than the universal asymptotic 
value, v = 0.588, due to the finite length of the chains. 




024680 2 4 



t k^ t k^ 

FIG. 2; Dynamic structure factor of the charged chain scaled 
with Zimm (left) and Rouse (right) exponent z for wavevec- 
tors in the range 0.2 < ka < 0.3 (large wavelengths) and 
times t < t^i^^. The inset shows the corresponding data for 
the uncharged chain. 



The inset of Fig. [T] shows the mean square displace- 
ment of the central monomer in the reference frame of 
the chain, ((Ar)^) = ((r/t)-R_(t)-r,(0)+R JO))^), as 
a function of time for charged and uncharged chains. If 
the time is scaled with the diffusion constant D, the two 
curves lie almost on top of each other, indicating that in- 
ternal modes are slowed down by the relaxation effect in 
the same way than the overall chain diffusion. The quan- 
tity ((Ar)^) reaches a plateau at the Zimm time tzi,„„^. 
From Fig. [Tl we infer Dt^imm^^cr^, in rough agreement 
with the heuristic estimate GDt^^^^ ~ R^—. This yields 
'•zimm'^^OOr for the charged chain, and t2i^,„^200T for 
the uncharged chain. 

To assess the dynamical properties of the chains, we 
study the dynamic structure factor S{k,t). In the dy- 
namic scaling regimeiiSt, it should exhibit the universal 
scaling behavior S{k,t) — S{k,0)f{tk'^) with the expo- 
nent z = 3 for Zimm chains and z — 2 + 1/1/ for Rouse 
chains. Hence the curves S{k,t)/ S{k,0) for different k 
vectors should collapse onto a single master curve f{tk^) 
when plotted against tk^ with the correct exponent z. 
Fig. [5] shows such scaling plots for large wavelengths, 
1/Rg < k < 0.3a^^, using f = 0.67. Not surprisingly, 
uncharged chains exhibit Zimm scaling. For charged 
chains, however, the data collapse much better in the 
Rouse scaling plot than in the Zimm scaling plot: Below 
the Zimm time, the chain behaves as if hydrodynamic 
interactions were absent. 

A closer inspection of Fig. [2] shows that this is not yet 
the full story. At very early times, the data collapse in 
the Zimm plot seems more convincing than in the Rouse 
plot. Indeed, we can identify a second characteristic time 
t„ ^ 55r, below which the data show Zimm scaling for 
a range of /c- vectors that includes both long and shorter 
length scales, l/R^ <k < 0.5a^^ (Fig.E]). For t>tg, the 
Zimm scaling breaks down for all k. For k < k^ ^ 0.3cr~^, 
one has a crossover to Rouse scaling. For k > k^, the data 
do not scale at all. 

To summarize, we find that hydrodynamic inter ac- 
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FIG. 3: Dynamic structure factor of the charged chain scaled 
with Zimm {z — 3, light circles) and Rouse {z = 3.5, black 
circles) for wavevectors 0.2 < ka < 0.5 and times t < tg = 55r 
(left) and t > tg (right). The inset on the right shows data for 
large wavelengths, k < kg = 0.3 /a at tg < t < t^.^^^ = 300r, 
which exhibit Rouse scaling as already demonstrated in Fig. 
[21 The main graph shows the data for short wavelengths, 
k > kg = 0.3/a and times to < t < 150r. 



tions are screened on large length scales. However, the 
screening is 'delayed' and not yet effective at early times 
t < tg 55r. A similar effect has been observed by 
Ahlrichs et al. in semidilute polymer meltsiS. In this 
case, the initial decay oiS{k,t) was found to be governed 
by an unscreened diffusion tensor. Screening sets in as 
soon as chains interact with each other, i.e., mesh blobs 
have moved their own size. Beyond that time, chain parts 
cannot follow the flow and act as immobile obstacles that 
produce Darcy-type friction. 

In our case, the phenomenon is similar, but the un- 
derlying physics is clearly different. The large reduc- 
tion of the diffusion constant in the presence of charges 
(Z?c < Dy) suggests that the chain dynamics is mainly 
driven by electrostatics, like in electrophoresis, where hy- 
drodynamic screening is also observed. However, the 
derivation of the screened Oseen tensor in electrophoresis 
relies crucially on the fact that the monomer and the sur- 
rounding ion cloud are subject to a homogeneous electri- 
cal field. Intrachain interactions produce inhomogeneous 
fields, thus the theory cannot be applied. Moreover, we 
need to explain the screening delay, i.e., the Zimm- Rouse 
crossover observed in our simulations. 

To analyze the problem in more detail, we begin with 
investigating the time scales governing the dynamics of 
charge distributions in our system. To this end, we cal- 
culate the dynamic "charge structure factor" Sq{k,t) = 
(|/g(k,t)/,(-k,0)|) with Iq = J2^<l^e'^'^^, where the 
sum i runs over all charged particles and q. denotes their 
charge. The result is shown in Fig. 2] (left). In the sim- 
plest linearized dynamical mean-field theory, the charge 
distribution decays exponentially towards the Debye- 
Hiickel distribution with the characteristic decay time 
T^{k) — l/(kgT{k'^ + A^-^). In the simulations, the de- 
cay of Sq{k,t) is slightly slower and depends less on fc - 
presumably due to the fact that the motion of the ions is 



FIG. 4: Illustration of time scales governing the charged sys- 
tem. Left: dynamic charge structure factor in the range 
kX^ < 1. Dashed lines show theoretically predicted expo- 
nential decay exp(— t/r(fc)) (see text) for A; = (upper curve) 
and k = (lower curve). Right: Autocorrelation function 

of Gi{t) and G2{t) (see text for definitions). Here ■ denotes 
dot product and : double-dot product (double contraction). 

not yet diffusive on the time scale . Nevertheless, the 
charge distribution has clearly relaxed on the time scale 
tg of the Zimm-Rouse crossover. On this time scale, the 
ions follow the chain adiabatically. 

We thus proceed to calculate the velocity field v gen- 
erated by the chain in the Debye-Hiickel and Stokes ap- 
proximation. For simplicity, we neglect the finite diam- 
eter of monomers and ions. The chain is characterized 
by the positions of the monomers a and their central 
charge . Each monomer charge generates a charge dis- 
tribution p{r), which in turn generates an electrostatic 
field E(r) and imparts a force field f.j (r) = p(r)E(r) to 
the fiuid. In addition, non-electrostatic interactions {e.g., 
bond interactions) between monomers a and (3 generate 
forces ('"0/3)^0/3 at r = R^, which are also transmitted 
to the fluid. Here r^^ = R^ — R^ and f = r/r. 

We consider the limit fc ^ l/A^, i.e., we focus on the 
length scale regime well above the Debye-Hiickel screen- 
ing length. After working out the theory, we obtain the 
velocity field v(k) = J2f3 exp(ik ■ R^) v('5)(k), 

with g{x) = e~^ (a;^/2 — x — 1) /x'^ and k — k/k, where kk 
denotes the tensor product. A derivation of this expres- 
sion shall be given elsewhere. From Eq. ([T]), one readily 
calculates the flow velocity at the position of a monomer 
a, v(R^) = X)/3 '^''^H'^/aa)- The contributions v^''^ are 
conveniently expanded in multipoles. Here we quote only 
the monopole and dipolc term 

vS'W^5;^i:{/.. + ^»(^)}(i + -)-.(2) 

The non-electrostatic interactions /^^ between 
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monomers only contribute to the monopole term, 
Eq. Inspecting the latter, one notices that it 

vanishes for locally straight chains, because contribu- 
tions of vectors r„ are then balanced by contributions 
of vectors r^^, — — r^^ ■ For symmetry reasons, the 
time-averaged conformations of the chain have this 
property. Therefore, one might speculate that the 
monopole becomes small compared to the higher order 
multipoles on sufficiently large time scales. On the 
other hand, it is clear from Eq. ([1]) that all multipoles 
must vanish on time scales where the vectors r„ are 
distributed isotropically. We thus need to compare the 
characteristic decay times for the monopole and the 
higher order multipoles quantitatively. To this end, we 
have evaluated the autocorrelation function of the vector 
= J^-ylfil-idiri^^/Xoyis^, which characterizes the 
monopole term, and that of the traceless part G^it) of 
the tensor J^^ <li3l-,9{'r I ^^d)^ ,}-,r ^^^^ which characterizes 
the dipole term. Here (3 was chosen to be the central 
monomer in the chain. The results are shown in Fig. 0] 
(right). For Gi, the decay is much faster than that for 
G2. By the time to(55r), the function {bfG^{t) ■ Gi(0)} 
has dropped by two orders of magnitude. This suggests 
that the dynamics of the chain is indeed governed by the 
higher order multipole interactions at later times. In the 
Stokes approximation, the latter vanish too, but much 
slower than the multipole term. If one includes nonlinear 
hydrodynamics, additional interactions come into play 
that presumably do not vanish at all. For example, the 
inhomogeneous velocity field created by the monomer 
/3 in the vicinity of monomer a convects the ion cloud 
surrounding a. This contributes an effective additional 
force field in the vicinity of , which in turn retroacts 
on /3, leading to an induced hydrodynamic interaction 
that remains finite at all times. 

Based on these considerations, we rationalize our sim- 
ulation results as follows: On short time scales, the 
monopole term, Eq. ([2]) dominates, hydrodynamic inter- 



actions are effective and lead to Zimm scaling. On longer 
time scales, time-averaged chains become locally straight 
and higher order terms take over. Zimm scaling is then 
replaced by Rouse scaling^ ^. The relevant characteris- 
tic length scale for the process of local chain straight- 
ening is twice the maximal distance between interacting 
monomers - in our case, roughly A,, ^ 3.4cr, since the 
mean distance between charged monomers in the chain 
is roughly qq — 1.7 a and A^^ < aq. This corresponds 
well with the observation that the Rouse scaling breaks 
down for large wave vectors with k < kQ 1/ Xq. It also 
explains why the crossover time to from Zimm to Rouse 
behavior does not seem to depend on the wavevektor k. 

To summarize, we have studied the dynamical behavior 
of polyelectrolytes in salt solution on time scales below 
the Zimm time. At very short times, they behave like 
regular self-avoiding chains that are subject to hydrody- 
namic interactions. In an intermediate time scale range 
between a microscopic crossover time and the Zimm time, 
however, the chain behaves like a Rouse chain, i.e., hy- 
drodynamic interactions are effectively screened. These 
findings should have implications for the interpretation of 
various processes of biological and nanotechnological in- 
terest. Many dynamical processes where polyelectrolytes 
are involved, such as DNA translocation through pores 
or RNA folding etc. are driven by internal chain modes. 
Our results indicate that the influence of hydrodynamics 
on these processes might be much smaller than is com- 
monly assumed. Theoretical models that neglect hydro- 
dynamics may be closer to reality than more refined mod- 
els that include hydrodynamics. 
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